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Abstract 

We report on evidence for black hole spindown in the light curves of the BATSE catalogue of 1491 long GRBs by 
application of matched filtering. This observation points to a strong interaction of the black hole with surrounding high 
density matter at the ISCO, inducing non-axisymmetric instabilities sustained by cooling in gravitational wave emission. 
Opportunities for LIGO- Virgo and the recently funded KAGRA experiments are highlighted, for long GRBs with and 
without supernovae and for hyper-energetic core-collapse supernovae within a distance of about 35 Mpc in the Local 
Universe. 
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1 Introduction 

Gamma-ray bursts (GRBs) and core-collapse super- 
novae (CC-SNe) are by far the most energetic and 
enigmatic transient events associated with neutron 
stars and mass black holes. A key objective is to 
identify the nature of their explosion mechanism. 

Extremely powerful events are unlikely powered 
by their presumably prodigious output in MeV neu- 
trinos. They may, instead, be powered by rota- 
tion (e.g. g). GRB 030329/SN 2003dh and GRB 
031203/SN 20031w are hyper-energetic events, requir- 
ing an energy reservoir that exceeds the maximal ro- 
tational energy of a proto-neutron star (PNS) by an 
order of magnitude [22] ■ These anomalous super- 
novae - of "Type I X" - practically rule out PNS as 
universal inner engine to all core-collapse supernovae. 

In light of these hyper-energetic events and the di- 
versity in GRBs in durations and associations (with 
and without supernovae, e.g., listed in [27]), we here 
turn to inner engines hosting a rotating black hole 
surrounded by high-density matter. These systems 
appear naturally in mergers of neutron stars with an- 
other neutron star or companion black hole, as well 
as in core-collapse of relatively mass stars. 

1.1 Frame dragging induced outflows 

According to the Kerr metric, the angular momen- 
tum of rotating black holes induces frame dragging. 
Frame dragging has recently been measured by the 
LAGEOS II and Gravity Probe B [TU] satellite ex- 
periments around the Earth, manifest in an angular 
velocity uj ~ J/r 3 at an orbital radius r and given 



the Earth's angular momentum J. By this scaling 
property, these experiments equivalently measured u 
around a maximally spinning Kerr black hole at a 
distance of about 5.3 million Schwarzschild radii. 

As part of the gravitational field, the frame drag- 
ging angular velocity u is a clean and universal agent 
inducing non-thermal radiation processes in space- 
time around rotating black holes. Relativistic frame 
dragging is encountered as uj approaches the angular 
velocity SIr of the black hole itself. 

Frame dragging around rapidly rotating black 
holes enables transfer of a major fraction of its spin 
energy to surrounding matter via an inner torus mag- 
netosphere simultaneously with transfer of a minor 
fraction to infinity in the form of ultra-relativistic 
outflows along the spin axis. The first can be seen by 
an equivalence in poloidal cross-section of the inner 
and outer faces of a torus in suspended accretion to 
the magnetospheres of neutron stars jTHl The 
second is described by a potential energy £ = ujJ p 
for test particles with angular momentum J p . Here, 
£ assumes energies on the order of Ultra High Energy 
Cosmic Rays (UHECRs) for particles in superstrong 
magnetic fields, typical for the catastrophic events 
under consideration [25j . Frame dragging is hereby a 
causal mechanism for the onset of a two-component 
outflow, comprising a baryon-rich wind from an inner 
disk or torus collimating the latter in the form of a 
baryon-poor jet. 

A two-component outflow with different baryon 
loading is advantageous for GRB-supernovae, whose 
GRB emissions are generally attributed to dissipa- 
tion of ultra-relativistic baryon-poor outflows and 
whose supernovae can be attributed to irradiation 
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from within [23) . such as by impulsive momentum 
transfer of internal magnetic winds onto the remnant 
stellar envelope. The efficiency of the latter favors 
relatively baryon-rich winds [271 , which generally will 
produce asphcrical explosions. These two-component 
outflows circumvent the limitations of neutron stars 
with approximately uniform baryon-loading through- 
out their wind, rendering these less amenable to 
making hyper-energetic supernovae with a successful 
GRB. 

1.2 If not a pulsar, how to identify 
rotation? 

In our model, the Tgo durations of long GRBs is at- 
tributed to the lifetime of rapid spin of the black 
hole. It sets the duration for energy and angular 
momentum transfer from the black hole mostly onto 
the surrounding matter. When the black hole has 
slowed down sufficiently, this interaction ceases, and 
the inner accretion disk or torus will fall in, heralding 
the end of the GRB. That is, a long GRB ends with 
hyper-accretion onto a slowly spinning black hole, 
whose angular velocity is approximately equal to that 
of the ISCO as prescribed by the Kerr metric. 

Furthermore, our model considers the aforemen- 
tioned two-component outflow. Any time-variability 
in the collimated torus wind - in strength or orien- 
tation - inevitably modulates the inner jet producing 
the GRB with possibly quasi-periodic variations in 
the light curve. The frequency scale of these vari- 
ations are those associated with the ISCO around 
the black hole, i.e., on the order of 1 kHz, possibly 
higher if multipole mass moments are involved and 
possibly lower by precession. A search for such high 
frequency QPOs requires the highest sampling fre- 
quencies available for GRB light curves, well beyond 
the 64 ms sampling interval in the BATSE catalogue. 

Following the above, we consider the problem of 
tracking the evolution of black hole spin, as deter- 
mined by its interaction with a surrounding torus. 
Attributed to frame dragging as indicated above, the 
count rate of gamma-ray photons observed should 
show a near-exponential decay as the angular veloc- 
ity approaches the fixed point fin = Qisco- 

The evolution of an initially rapidly rotating black 
hole that ensues is mostly in a decrease in its angu- 
lar momentum and less so in a decrease of its total 
mass [2U [26]. Since oj and hence £ decrease along 
with fin, the luminosity of the baryon-poor outflows 
will decrease in the process of black hole spindown in 
approaching a fixed point fin = Qd- Here, Flo de- 
notes the angular velocity of the inner disk, which is 
expected to closely track the angular velocity Qisco 
of the Inner Most Stable Circular Orbit (ISCO), as 
described by the Kerr metric. Thus, an asymptoti- 
cally exponential decay in the light curve of any high 



energy emission produced by dissipation of the afore- 
mentioned baryon poor outflows is a characteristic 
imprint of black hole spindown. 

The above suggests extracting a normalized light 
curve (nLC) of long GRBs, to identify the anticipated 
late time, exponential decay. Here, we report on this 
study using a recent implemented matched filtering 
[25] , to identify a model light curve in the complete 
BATSE catalogue of 1491 long GRBs. Matched fil- 
tering allows for an accurate extraction of a normal- 
ized light curve, in which fluctuations on all interme- 
diate time scales are filtered out, that can be used 
to validate model templates. We apply this to vari- 
ous model templates representing spindown, of black 
holes and (proto-)neutron stars. 

1.3 EM priors to GW searches 

Properties of the inner engine of hyper-energetic su- 
pernovae and long GRBs in the electromagnetic spec- 
trum of radiation can provide valuable priors to 
searches for their gravitational wave emissions. By 
the very nature of these catastrophic events, the high 
density matter at nuclear densities orbiting at about 
the Schwarzschild radius of the central engine, is ex- 
pected to develop a non-axisymmetric mass distribu- 
tion in response to exceptionally large energy fluxes. 
If so, significant emission in gravitational waves is in- 
evitable for a sustained period of time, provided by 
a balance between heating, apparent in MeV neu- 
trino emission, and cooling, by aforementioned mag- 
netic disk winds and gravitational radiation. This 
outlook offers some novel opportunities for the ad- 
vanced detectors LIGO- Virgo and KAGRA currently 
under construction, that are expected to be opera- 
tional within this decade. 



normal 




mean =4.72e+003 
a = 659 









3000 3500 4000 4500 5000 5500 6000 6500 7000 



v [km/s] 



1 

- mean =5.68e+003 
o = 860 

















3000 3500 4000 4500 5000 5500 6000 6500 7000 



v [km/s] 

Figure 1: Shown are the histograms of the ejection 
velocities of the sample of [14] of CC-SNe with nar- 
row (top) and broad (bottom) line emission lines. The 
mean of the latter is larger than that of the former 
by 3.9cr. 

Of particular interest are priors to identifying 
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black holes or PNS, in view of their dramatically dif- 
ferent outlook in gravitational wave emission. The 
latter offers a broad range of possible radiation chan- 
nels which, however, remain hitherto somewhat un- 
certain. Examples are acoustic modes, convection, 
differential rotation, and magnetic fields [17j HH El [8] 
[H [TTJ [SJ [T] [15] . If an electromagnetic prior is found 
to rule out a PNS, we can direct our attention to 
black holes as a probable alternative. An attractive 
possibility, for example, is a prolonged quasi-period 
emission by quadrupole mass inhomogeneities or non- 
axisymmetric instabilities in the inner disk or torus 
at twice the orbital frequency, particularly so if pow- 
ered by the spin energy of the black hole. 

Table 1: Selected supernovae with kinetic energies 
Esn in units of 10 51 erg. Required energy reservoir 
E rot is expressed in E = E rot /E c , E c = 3 x 10 52 erg 
(adapted from [27]). 
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2 Diversity out of universality 

In contrast to Type la supernovae, CC-SNe are not of 
one kind with both narrow and broad emission lines, 
shown in the histograms of ejection velocities (Fig. 
1) of a sample of CC-SNe compiled by [14]. Those 
of broad emission line events are higher on average 
that those of narrow emission line emission events 
with a statistical significance in excess of 4 a. The 
explosions of the former appear to be relatively more 
energetic. 

Rapidly rotating black holes (Qh > QiscOi 
a/M > 0.36) can form in CC-SNe and mergers of 
two neutron stars. If the progenitor of the former is 
a member of a short period binary, a rapidly rotating 
black hole is formed with a rotational energy between 
about one- and two-thirds of the extremal value for 
Kerr black holes. In particular, dimensionless specific 
angular momenta a/M = 0.7679 and a/M = 0.9541 
are found with E rot /E™t x — 0.3554 and, respec- 
tively, E rot /E™t x = 0.6624. Neutron star- neutron 
star mergers form rapidly rotating black holes of rela- 
tively small mass, M, essentially the sum of the mass 
Mats of individual neutron stars except for a minor 



loss of mass that forms an accretion disk, for which 
we estimate 

2 / R 

— = - J — = 0.5963, 0.7303, 0.8433 (1) 

M 5 V R g ' ' w 

for Mns = 3, 2, 1.5M©. These values are in remark- 
able agreement with numerical simulations, showing 
a/M = 0.74 - 0.84 for M NS = 1.5M Q [3]. In the 
merger of a neutron star with black hole compan- 
ion, an accretion disk will form from tidal break-up 
around the latter if the black hole is not too massive 
(the limit for which is relaxed if it spins rapidly). In 
this event, the spin of the black hole is unchanged 
from its spin prior to the merger. Consequently, the 
rotation of the black hole may be diverse, given by 
the diversity spin in neutron star-black hole binaries. 

In considering GRBs from rotating black holes, 
we associate long/short GRBs with rapidly/slowly 
spinning black holes. The above leads to the perhaps 
counter-intuitive conclusion that mergers of neutron 
star-neutron star binaries produce long GRBs, more 
likely so than short GRBs [24] . They may further 
be produced by mergers of neutron stars with black 
hole companions (e.g. [UlIU]), especially those with 
rapid spin (T9"] . 

Consequently, long GRBs are expected to occur 
both with and without supernovae involving rapidly 
rotating black holes. Notable examples are given in 
Table 1 of 127], e.g., GRB 060614 with T 90 = 102 s. 

3 A light curve of long GRBs 

To leading order, the evolution of rotating black hole 
interacting with surrounding matter via a torus mag- 
netosphere is described by conservation of energy and 
angular momentum associated with the black hole 
luminosity L#, its mass M and angular momentum 
Jh , satisfying [19] 

L H = -M, J H = -T, (2) 

where Lh = QrT associated with the torque T = 
k(Qh — Qd)- Here, k incorporates the physical and 
geometrical properties of the inner torus magneto- 
sphere, and is taken to be tightly correlated to 
&ISCO mentioned earlier. Upon numerical integra- 
tion of ([2]), a model light curve can be calculated 
from the resulting baryon-poor outflow as a function 
of Qh = ^if(i) and 9h = 0ff(t) in its dependence on 
the radius Rd = Risco of the inner disk. For our 
intended leading order analysis, we assume that the 
photon count rate is proportional to the luminosity in 
the baryon poor outflow. Fig. 2 shows the resulting 
template, scaled to a duration of about one minute. 

A finite polar cap supports an open magnetic flux 
tube to infinity, while the remainder of the black 
hole event horizon is connected to the inner disk or 
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torus via an inner torus magnetosphere. The time- 
averaged luminosity of the jet along the open mag- 
netic flux tube that results from the action of frame 
dragging, satisfies the angular scaling [23] 

with further dependence on the angular velocity of 
the black hole according to Lj oc Vl H z 2 Ek,T , Ek,D °c 

(n D R D ) 2 e(z), e(z) = Jl-£, Rt = zM. Here, we 
consider the geometrical ansatz O 4 ^ oc z n with n = 2. 
(Here, n = 2 corresponds to the expansion of the sur- 
face area Ah — A\z + A2Z 2 + • • • of the polar cap as a 
function of the radius of the inner disk.) Fig. 2 shows 
the model light curve produced by ([3]) following ([2]). 




a/M (e H /max(9 H ) 



Figure 2: The template light curve in gamma-rays 
for an initially extremal black hole, shown as a func- 
tion of time (top left), dimensionless angular mo- 
mentum (bottom left), and the associated horizon 
half-opening angle Oh (right). During initial spin- 
down, the rise is due to a marked increase of Oh with 
the expansion of the ISCO, leading to a maximal lu- 
minosity at about 16% of the duration of the burst. 
There is an accompanying decrease in the total mass 
of the black hole [24l|27]. 

Thus, the light curve ((2J is based on two posi- 
tive correlations, in the GRB luminosity with uj and 
Oh with the radius of the ISCO. For the latter we 
use the Ansatz n=2, a choice among a few integers 
by insisting on analyticity. We now seek validation 
of this Model A in the BATSE light curves of long 
GRBs by application of matched filtering. 

4 Normalized GRB light curve 

To extract an intrinsic light curve of long GRBs, 
we set out to filter out all fluctuations, randomly or 



quasi-periodic, in the BATSE data (Fig. 3) that are 
not representative for the evolution of the black hole. 
These include, without being exhaustive, precession 
of the disk or black hole and instabilities, in the ac- 
cretion disk, the inner torus magnetosphere and in 
the interface between baryon-poor jet and the sur- 
rounding baryon-rich disk winds, as well as fluctua- 
tions in the dissipative fronts due to turbulence. To 
this end, a normalized light curve (nLC) is extracted 
by an ensemble average of individually normalized 
light curves defined by a best-fit to a template upon 
scaling and translation of time and count rate. 
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Figure 3: Compilation of the complete BATSE cata- 
logue of 1491 light curves of long GRBs, sorted by 2 
s < Tgo < 1307 s. Each light curve shown represents 
the sum of the photon count rate in all four BATSE 
energy channels, is smoothed with a time scale of 2.56 
s and is plotted as a function of time normalized to 
T 90 . (Reprinted from [28].) 



Fig. 4 shows the results along with the residuals, 
showing convergence for very long bursts (Tqq > 20 s) 
to Template A (Fig. 2). The results are sub-optimal 
for other model templates, here Template B for spin- 
down against surrounding matter beyond the ISCO 
satisfying Sl^i = ^h, and Template C, for spindown 
of a PNS. 

The preferred match to Template A establishes 
confidence in the natural in view of spindown against 
matter at the ISCO, that represents a suspended ac- 
cretion state with a luminous output in gravitational 
radiation [20]. A detailed analysis points to a sen- 
sitivity distance of about 35 Mpc by the advanced 
gravitational wave detectors LIGO- Virgo and KA- 
GRA [20], by extrapolating the sensitivity distance 
determined from model injections into the strain am- 
plitude noise of the TAMA 300, shown in Fig. 5. 
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Figure 4: Shown are the nLC (thick lines) generated by model templates AC (thin lines) for the ensemble 
of 531 long duration bursts with 2 s < Tg < 20 s (left) and the ensemble of 960 long bursts with T 90 >20 s 
(right) and the associated deviations for Templates AC. Here, the standard deviation a is calculated from the 
square root of the variance of the photon count rates in the ensemble of individually normalized light curves as 
a function of normalized time. (Reprinted from 28 .) 
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Figure 5: The long duration negative chirp can be 
searched for using a dedicated time-sliced matched 
filtering procedure, here illustrated following signal 
injection in TAMA 300 noise data. It points to a 
sensitivity distance of about 35 Mpc for the upcom- 
ing advanced detectors. (Reprinted from |26|.) 



5 Discussion and conclusions 

Some hyper-energetic supernovae - of Type IX- 
appear to require an energy reservoir that exceeds 
the maximal rotational energy of a PNS. GRB- 
supernovae ideally has an inner engine that pro- 
duces simultaneously a baryon-rich wind to power 
the supernova and a baryon-poor jet to power the 
GRB. These and other considerations point to a uni- 
versal inner engine in the form of a rotating black 
hole, whose diversity in spin [19] can account for 
short/long GRBs from hyper- and suspended accre- 
tion [2"3"] . 

The Fanaroff-Riley class I and II radio-loud AGN 
may represent a closely related dichotomy for super- 
massive black holes, where the latter may be related 



to a similar suspended accretion state catalyzing 
black hole spin energy into powerful winds. eLISO 
may test for suspended accretion by searches for 
an accompanying low-frequency emission in gravita- 
tional waves, notably from SgrA*. The Type B X-ray 
flaring in the galactic microquasar GRS 1915+105 
may further illustrate the suspended accretion state 
in a galactic microquasar |22j . 

Fig. @] validates our model for black hole spin- 
down most likely against matter at the ISCO, rather 
than further out or the spindown of (proto-)neutron 
stars. As a corollary, some of the neutron star- 
neutron star mergers are expected to long GRBs, 
more likely so than short GRBs. 

Our results point to specific high-frequency chirps 
both in GWs of interest to LIGO- Virgo and KAGRA, 
and in intensity modulations in prompt GRB emis- 
sions to some of the gamma-ray satellite missions. 
For the first, we estimate a sensitivity distance to be 
about 35 Mpc for the advanced generation of gravi- 
tational wave detectors, indicated in Fig. [5l 
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DISCUSSION 

J. Beal: Is there anisotropy in the emitted gravita- 
tional radiation? 

M. van Putten: The anisotropy in quadrupole 
emissions is known to be small, about a factor of 
1.58 at most. 

B. Czerny: (a) GW efficiency will depend strongly 
on the m in azimuthal modes. Are low m more im- 
portant in the Papaloizou-Pringle instability? (b) 
Why don't you use the term "hyper nova" ? (c) How 
sensitive are your GRB templates to the detailed as- 
sumptions that go into BZ efficiency? 

M. van Putten: (a) In [5], we considered the prob- 
lem of the spectrum in gravitational wave emission 
using a toy model for a magnetized disk. The re- 
sults indicate a dominant emission in m = 2. In |21) . 
I extended the PP instability analysis to tori with 
finite width, and found that the low m modes (start- 
ing with m = 1) become unstable first in response 
to thermal or magnetic pressures. Both results are 
fortunate, pointing to major emission in quadrupole 
radiation, that falls within the limited band width of 
sensitivity of the LIGO- Virgo and KAGRA detectors. 

(b) In seeking EM priors on the inner engine of 
CC-SNe, I consider the possibility of ruling out PNS 
by determining the required energy reservoir relative 
to E c = 3 x 10 52 erg, the canonical value for the 
maximal energy in rotation of a PNS. The term "hy- 
pernova" has been introduced for supernovae with 
isotropically equivalent kinetic energies of a few times 
10 52 erg, well in excess of typical SN energies of a few 
times 10 51 erg. This definition is not sufficiently spe- 
cific to serve as a prior on the inner engine. 

(c) BZ proposed an open model for a black hole 
luminosity entirely in open outflows, sustained by ac- 
cretion of magnetized flows with no feedback onto 
the inner accretion disk. The black hole evolves ac- 
cording to the net result of losses (in energy and an- 
gular momentum) in magnetic outflow and gain by 
accretion. Hyper-accretion will likely cause the black 
hole to spin up continuously up to close to maximal 
rotation [13] . In associating the duration of hyper- 
accretion with the observed Tqq of the bursts, a model 
light curve from a baryon poor jet along the black 
hole spin axis should be increasing, followed by a rel- 
atively sharp drop off when accretion ceases. This is 
at odds with the gradual and extended decay seen in 
the nLC obtained by Templates A-C (Fig. 4). 



